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The masses of the heavy tetraquarks with open charm and bottom are calculated 
within the diquark-antidiquark picture in the framework of the relativistic quark 
model. The dynamics of the light quarks and diquarks is treated completely relativis- 
tically. The diquark structure is taken into account by calculating the diquark-gluon 
form factor. New experimental data on charmed and charmed-strange mesons are 
discussed. Our results indicate that the anomalous scalar L>* (2317) and axial vector 
Z) s i(2460) mesons could not be considered as diquark-antidiquark bound states. On 
the other hand, _D S (2632) and D*j(2860) could be interpreted as scalar and tensor 
tetraquarks, respectively. The predictions for masses of the corresponding bottom 
r~| \ counterparts of the charmed tetraquarks are given. 



PACS numbers: 12.40.Yx, 14.40. Gx, 12.39.Ki 



INTRODUCTION 



In last few years a significant experimental progress has been achieved in charmonium 
and charmed meson spectroscopy. Many new states, such as X(3872), y(4260), Y(4360), 
y(4660b Z(4248), Z(4430), £>* (2317), L> s i(2460), £> s (2632) and £)* </ (2860) etc., were ob- 
served pj which cannot be simply accommodated in the quark-antiquark (cc and cs) picture. 
These states and especially the charged charmonium-like ones can be considered as indica- 
rS \ tions of the possible existence of exotic multiquark states. 

The open charm mesons, both with and without open strangeness, represent a special in- 
terest. Even seven years after the discovery of the charmed-strange _D* (2317) and -D s i(2460) 
mesons their nature remains controversial in the literature. The abnormally light masses of 
these mesons put them below DK and D*K thresholds, thus making these states narrow, 
since the only allowed strong decays D^J — > D^tt violate isospin symmetry. The pecu- 
liar feature of these mesons is that they have masses almost equal or even lower than the 
masses of their charmed counterparts -Dq(2400) and -Di(2427) [2j-y|. Most of the theoretical 
approaches including lattice QCD 0, QCD sum rule p] and different quark model 0, [1] 
calculations give masses of the + and 1 + P-wave cs states significantly heavier (by 100-200 
MeV) than the measured ones. Different theoretical solutions of this problem were proposed 
including consideration of these mesons as chiral partners of 0~ and 1~ states [3], cs states 
which are strongly influenced by the nearby DK thresholds DK or D s tt molecules (llj . 
a mixture of cs and tetraquark states T3"3]. However the universal understanding of their 
nature is still missing. Therefore it is very important to observe their bottom counterparts. 



The unquenched lattice calculations of their masses can be found in Ref. [17 . 
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TABLE I: Predictions 24| for the masses of charmed 
MeV) in comparison with available experimental data 
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new states recently observed by BaBar [25J] 



Another unexpectedly narrow charmed-strange meson £) s (2632) was discovered by SE- 
LEX Collaboration Il8( . Its unusual decay properties triggered speculations about its pos- 
sible exotic orig in |l|. However, the status of this state remains controversial since FOCUS 
and Belle ED reported negative results in their search for this state. 



19 , BaBar 20 



Three other charmed-strange mesons (2710), _D*j(2860) and _D s j(3040) were discov- 



ered at 5-factories by Belle and BaBar |22|, |23J. Not all of them could be simply accom- 



modated in the usual cs picture. Their decay pattern implies that _D*j(2860) should have 
natural parity, while D s i(2710) and _D s j(3040) should have unnatural parity. Our recent cal- 
culation of the heavy-light quark- ant iquark meson spectra [24J have shown that £) sl (2710) 
and -D s j (3040) are good candidates for the 2 3 5*i and 2P 1 states, respectively. They nicely 
fit to the corresponding Regge trajectories, while -D* (2317), -D s i(2460), .D*j (2860) and 
_D S (2632) have anomalously low masses (lower than expected by 100-200 MeV) and do not 
lie on the respective Regge trajectories. 

In the charmed sector very recently the BaBar Collaboration discovered four new signal 
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peaks £>(2550), L>*(2600), L>(2750) and £>*(2760) |25|. The last two signals, observed in 
D*7r and Dtt modes, have close mass and width values (which differ by 2.6a and 1.5a, 
respectively) and therefore could belong to the same state. The angular analysis shows that 
these signals can be considered as candidates for the radially excited 2 1 S , 2 3 Si states and 
orbitally excited 1 3 D\ state, respectively. Their mass values are in a good agreement with 
the results of our model. We summarize our predictions for the masses of the cq (q = u, d) 
and cs mesons in Table Hand confront them with the available experimental data [4], |25[. As 
it is clearly seen there are indications of the possible existence of the exotic states especially 
in the charmed-strange sector. 



In our papers [261 . |27| we calculated masses of the hidden heavy-flavour tetraquarks and 
light tetraquarks in the framework of the relativistic quark model based on the quasipotential 
approach in quantum chromodynamics. Here we extend this analysis to the consideration 
of heavy tetraquark states with open charm and bottom. This study could help in revealing 
the nature of the anomalous charmed-strange mesons. 

As previously (26|, 27[, we use the diquark-antidiquark picture to reduce a complicated 
relativistic four-body problem to two subsequent more simple two-body problems. The first 
step consists in the calculation of the masses, wave functions and form factors of the diquarks, 
composed from light and heavy quarks. At the second step, a tetraquark is considered to be 
a bound diquark-antidiquark system. It is important to emphasize that we do not consider 
the diquark as a point particle but explicitly take into account its structure by calculating 
the form factor of the diquark-gluon interaction in terms of the diquark wave functions. 

II. RELATIVISTIC MODEL OF TETRAQUARKS 

In the quasipotential approach and diquark-antidiquark picture of heavy tetraquarks 
the interaction of two quarks in a diquark and the diquark-antidiquark interaction in a 
tetraquark are described by the diquark wave function (^d) of the bound quark-quark state 
and by the tetraquark wave function (\I/t) of the bound diquark-antidiquark state, respec- 
tively. These wave functions satisfy the quasipotential equation of the Schrodinger type 

H 

where the relativistic reduced mass is 

E l E 2 M 4 - (m 2 - m 2 ) 2 

m ~ e^+W 2 ~ up ' (2) 

and Ei, E 2 are given by 

_ M 2 -m\+m{ _ M 2 -m\ + m 2 2 

1= 2M ' 2= 2M ' ^ 

Here, M = Ei + E 2 is the bound-state mass (diquark or tetraquark), m 12 are the masses 
of quarks (q and Q) which form the diquark or of the diquark (d) and antiquark (cf ) which 
form the heavy tetraquark (T), and p is their relative momentum. In the center-of-mass 
system the relative momentum squared on mass shell reads 

2 [M 2 -(m l + m 2 ) 2 ][M 2 -(m 1 -m 2 ) 2 ] 

1 ' 4M 2 ' 1 ' 
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The kernel V(p, q; M) in Eq. ([T]) is the quasipotential operator of the quark-quark or 
diquark-antidiquark interaction. It is constructed with the help of the off-mass-shell scatter- 
ing amplitude, projected onto the positive-energy states. For the quark-quark interaction in 
a diquark we use the relation V qq = V qq -/2 arising under the assumption of an octet structure 
of the interaction from the difference in the qq and qq colour states. An important role in 
this construction is played by the Lorentz structure of the confining interaction. In our anal- 
ysis of mesons, while constructing the quasipotential of the quark-antiquark interaction, we 
assumed that the effective interaction is the sum of the usual one-gluon exchange term and a 
mixture of long-range vector and scalar linear confining potentials, where the vector confin- 
ing potential contains the Pauli term. We use the same conventions for the construction of 
the quark-quark and diquark-antidiquark interactions in the tetraquark. The quasipotential 
is then defined as follows 0, 26 



(a) For the quark-quark (Qq) interactions, V(p, q; M) reads 



V(p, q; M) = Ui(p)u 2 (-p)V(p, q; M)u 1 (q)u 2 (-q), 



(5) 



with 



V(p,q;M) 



a a £Uk)-rf 7 £ + Konf(k)r?(k)r 2;/i (-k) + Kf nf (k) 



Here, a s is the QCD coupling constant; Dn V is the gluon propagator in the Coulomb gauge, 



D 00 (k) = - 



47T 



D ij (k) 




k { k j 



D 



Oi 



D 



iO 



(6) 



and k = p — q; 7 M and u(p) are the Dirac matrices and spinors 

u\p) 



e(p) + m 



\ Hp) 



i 

cr ■ p 



X 



(7) 



V e (p) + m 



with e{p) = Vp 2 + m 2 - 

The effective long-range vector vertex of the quark is defined by 0] 



■IK 



r„(k) = 7m + ^<v r > 



fe = (0,k), 



(8) 



where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. In configuration space the vector and scalar confining potentials in the 
nonrelativistic limit |28| reduce to 



KonfW = (l-^)Konf(r), 
KonfW = ^conf(r), 

with 

Konf(r) = V c s on{ (r) + V£Jr 

where e is the mixing coefficient. 

(b) For the diquark-antidiquark (dd') interaction, V(p, q; M) is given by 



Ar + B, 



(9) 
(10) 



ZyEdEd 3 2y/Ed>Ed> 



+r d (PWAP') ^^Konf(k) + Konf(k) MQ)MQ') 



(ii) 
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where {d(P)\J fl \d(Q)) is the vertex of the diquark-gluon interaction which takes into account 
the finite size of the diquark [pO = {E d o , ±p) and = (E d( > h ±q), E d = (M 2 - Mj, + 

MJ)/(2M) and E d > = (M 2 - Mj + M|)/(2M) ' . 

The diquark state in the confining part of the diquark-antidiquark quasipotential fTTTj) is 
described by the wave functions 




for a scalar diquark, 

for an axial- vector diquark, 



(12) 



[Gd ■ PJP 



M d (E d (p) + M d 



(13) 



2 + M 2 



is the polarization vector of the axial- vector diquark with momentum p, E d (p) = y p -r j-v± d , 
and e d (0) = (0, e d ) is the polarization vector in the diquark rest frame. The effective long- 
range vector vertex of the diquark can be presented in the form 



J, 



2 



(/;/( 



E d E d 

(P + Q)u 



+ 



2JE d E d 



^pd yi/ 1, 



for a scalar diquark, 



for an axial-vector diquark. 



(14) 



Here, the antisymmetric tensor is defined by 



pa, 



(15) 

-iSkiu l^d is the total chromo- 



and the axial-vector diquark spin S d is given by (S d -k)u 
magnetic moment of the axial- vector diquark. 

The constituent quark masses m c = 1.55 GeV, = 4.88 GeV, m u = m d = 0.33 GeV, 
m s = 0.5 GeV and the parameters of the linear potential A = 0.18 GeV 2 and B = —0.3 GeV 
have values typical in quark models. The value of the mixing coefficient of vector and scalar 
confining potentials e = — 1 has been determined from the consideration of charmonium 
radiative decays [28| and the heavy-quark expansion [2^]. The universal Pauli interaction 
constant k — — 1 has been fixed from the analysis of the fine splitting of heavy quarkonia 
3 Pj - states [28]. In this case, the long-range chromomagnetic interaction of quarks vanishes 
in accordance with the flux-tube model. 

Since we deal with diquarks and tetraquarks containing light quarks and diquarks, re- 
spectively, we adopt for the QCD coupling constant a s (/i 2 ) the simplest model with freezing 
namely 



4:71 



A) In 



// 



2 + M 2 



A) 



n - - 



71/ , 



(16) 



B 



A 2 



where the scale is taken as [i = 2m 1 m 2 /(m 1 +m 2 ), the background mass is Mb = 2.24 
0.95 GeV S, and the parameter A = 413 MeV was fixed from fitting the p mass [31] 



that the other popular parametrization of a s with freezing [32j leads to close values. 



A = 
Note 
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TABLE II: Masses M and form factor parameters of diquarks. 
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III. DIQUARK AND TETRAQUARK MASSES 

At the first step, we calculate the masses and form factors of the heavy and light diquarks. 
Since the light quarks are highly relativistic a completely relativistic treatment of the light 
quark dynamics is required. To achieve this goal, we closely follow our consideration of 
diquarks in heavy baryons and adopt the same procedure to make the relativistic potential 



local by replacing ei^Qo) = Jm\ 2 + P 2 — >• £u,2 = (M 2 — rn^ x + m\ 2 )/2M. Solving numeri- 
cally the quasipotential equation ([!]) with the complete relativistic potential, which depends 



on the diquark mass in a complicated highly nonlinear way [33|, we get the diquark masses 
and wave functions. In order to determine the diquark interaction with the gluon field, which 
takes into account the diquark structure, we calculate the corresponding matrix element of 
the quark current between diquark states. Such calculation leads to the emergence of the 
form factor F(r) entering the vertex of the diquark-gluon interaction (33[. This form factor 
is expressed through the overlap integral of the diquark wave functions. Our estimates show 
that this form factor can be approximated with a high accuracy by the expression 



F(r) = l-e-^ 2 . (17) 

The values of the masses and parameters £ and £ for light and heavy scalar diquark [• • •] and 
axial vector diquark {• • •} ground states were calculated previously 0, 33[ and are given in 
Table HH 

At the second step, we calculate the masses of heavy tetraquarks considered as the bound 
states of a heavy-light diquark and light antidiquark. For the potential of the S'-wave 
((L 2 ) = 0) diquark-antidiquark interaction ffTTj) we get [27[ 



V(r) = Vboui(r) + V coni (r) + -L- jp [v Coul (r) + V^(r)] p - ^^(r) 
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+ - 
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Ag E,E 2 
4 MiM 2 



conf V 



Si-S s 



(18) 



where 



VboulM 



4 Fi(r)F 2 (r) 

- — CK S 

3 r 



is the Coulomb-like one-gluon exchange potential which takes into account the finite sizes 
of the diquark and antidiquark through corresponding form factors F 12 (r). Here, Si )2 are 
the spin operators of diquark and antidiquark. In the following we choose the total chromo- 
magnetic moment of the axial- vector diquark jj,d = 0. Such a choice appears to be natural, 
since the long-range chromomagnetic interaction of diquarks proportional to fid then also 
vanishes in accordance with the flux-tube model. 

The resulting quasipotential equation with the complete kernel (|18|) is solved numerically 
without any approximations. 

IV. RESULTS AND DISCUSSION 



Masses of the heavy tetraquark ground (IS) states calculated in the diquark-antidiquark 
picture are presented in Table HTT1 In this table we also give possible experimental candidates 
for charmed and charmed-strange tetraquarks. 

In the charmed meson sector the situation is rather complicated. Comparing results 
presented in Tables HI and II III we see that our model predicts very close masses for the scalar 
+ orbitally excited cq meson (l 3 Po) and the ground state (IS") tetraquark, composed from 
the scalar [cq] diquark and scalar [qq] antidiquark. The same is true for the masses of the 
1 + axial vector cq meson (lPi) and the tetraquark, composed from the axial vector {cq} 
diquark and scalar [qq] antidiquark. The calculated masses are consistent with the measured 
masses of the scalar Dq(2400) and axial vector Dx(2A30) mesons. The mixing of the cq and 
tetraquark states could be responsible for the observed difference in masses of the charged 
and neutral Dq(2400) mesons. 

In the charmed-strange sector the + and 1 + tetraquarks are predicted to have masses sig- 
nificantly (by 200-300 MeV) higher than experimentally measured masses of the _D* (2317) 
and _D s i(2460) mesons (cf. Tables H] IHIj) . This excludes the interpretation of these anoma- 
lously light D s mesons as the heavy diquark-antidiquark (tetraquark) states in our model. 
Instead, we find that the lightest scalar + tetraquark, composed from the scalar [cq] di- 
quark and scalar [sq] antidiquark, has a mass consistent with the controversial _D S (2632) 
observed by SELEX |l8|. The D*j(2860) meson, observed by BaBar Q both in DK and 
D*K modes, 1 has mass coinciding within experimental error bars with the prediction for the 
mass of the tensor 2 + tetraquark, composed form the axial vector {cq} diquark and axial 
vector {sq} antidiquark. The rather high value of its spin can explain the non-observation 
of this state by Belle in B decays. 

In Table [IV] we compare our results for the masses of the charmed-strange diquark- 



antidiquark bound states with the predictions of Refs. [13l-ll6l|. From this table we see that 



only the model of Ref. [13| gives masses of scalar and axial vector tetraquarks compatible 



with the observed masses of the _D* (2317) and _D s i(2460) mesons. All other models predict 



1 This state should therefore have natural parity and total spin J > 1. 
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TABLE III: Masses of charmed and bottom diquark-antidiquark ground (IS) states (in MeV) and 
possible experimental candidates |4]. S and A denote scalar and axial vector diquarks. 
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significantly higher mass values for these states. The main difference between these ap- 
proaches consists in the substantial distinctions in treating quark dynamics in tetraquarks. 
The authors of Ref. [l3[ use a phenomenological approach, determining diquark masses 
and parameters of hyperfme interactions between quarks from adjusting their predictions to 
experimental observables. Contrary we describe diquarks and tetraquarks dynamically as 
quark-quark and diquark-antidiquark bound systems and calculate their masses and form 
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TABLE IV: Comparison of theoretical predictions for the masses of cqsq tetraquarks (in MeV). 



this paper [13] [14] [y] [16]_ 

0+ 2619 2371 2840 2731 2616 

1+ 2678 2410 2841 

1+ 2723 2462 2880 2841 

0+ 2689 2424 2503 2699 

1+ 2757 2571 2748 

2+ 2863 2648 2983 2854 



factors in the model where all parameters were previously fixed from considerations of meson 



properties. Different dynamical approaches were applied in Refs. [14-16]|. The authors of 
Ref. [lij calculate diquark-antidiquark mass spectra in the quark model employing the QCD 
potential found by means of the AdS/QCD correspondence. Tetraquark masses are calcu- 
lated in Ref. 15[ in the nonrelativistic quark model including both the confining interaction 
and meson exchanges, while in Ref.(l6| the coupled-channel formalism is employed. In the 
latter two approaches the anomalous charmed-strange mesons could be only accommodated 
as a mixture of quark-antiquark and tetraquark states with a phenomenologically adjusted 
mixing interaction. Thus it seems to be unlikely that the _D* (2317) and _D sl (2460) could be 
pure diquark-antidiquark bound states. 



V. CONCLUSIONS 

We calculated the masses of heavy tetraquarks with open charm and bottom in the 
diquark-antidiquark picture using the dynamical approach based on the relativistic quark 
model. Both diquark and tetraquark masses were obtained by the numerical solution of the 
quasipotential wave equation with the corresponding relativistic potentials. The diquark 
structure was taken into account in terms of diquark wave functions. It is important to 
emphasize that, in our analysis, we did not introduce any free adjustable parameters but 
used their values fixed from our previous considerations of heavy and light hadron proper- 
ties. It was found that the _Dg(2400), £> s (2632) and D*j(2860) mesons could be tetraquark 
states with open charm, while the _D* (2317) and D s i(2460) mesons cannot be interpreted 
as diquark-antidiquark bound states. The masses of the bottom counterparts of charmed 
tetraquarks were calculated. It is important to search for them in order to help revealing 
the nature of controversial charmed and charmed-strange mesons. 
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